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Abstract 

This article addresses the boundary layer flow and heat transfer in third grade fluid over an unsteady permeable stretching 
sheet. The transverse magnetic and electric fields in the momentum equations are considered. Thermal boundary layer 
equation includes both viscous and Ohmic dissipations. The related nonlinear partial differential system is reduced first into 
ordinary differential system and then solved for the series solutions. The dependence of velocity and temperature profiles 
on the various parameters are shown and discussed by sketching graphs. Expressions of skin friction coefficient and local 
Nusselt number are calculated and analyzed. Numerical values of skin friction coefficient and Nusselt number are tabulated 
and examined. It is observed that both velocity and temperature increases in presence of electric field. Further the 
temperature is increased due to the radiation parameter. Thermal boundary layer thickness increases by increasing Eckert 
number. 
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Introduction 

There is a substantial interest of the recent researchers in the 
flows of non-Newtonian fluids. Such motivation in these fluids is 
mainly because of their use in the industrial and technological 
applications. Many materials like mud, pasta, personal care 
products, ice cream, paints, oils, cheese, asphalt etc. are non- 
Newtonian fluids. Most biological fluids with higher molecular 
weight components are also non-Newtonian in nature. The usual 
properties of polymer melts and solutions together with the 
desirable attributes of many polymeric solids, have given rise to the 
world-wide industry of polymer processing. The non-Newtonian 
fluids in particular have key importance in geophysics, chemical 
and nuclear industries, material processing, oil reservoir engineer- 
ing, bioengineering and many others. Rheological properties of all 
the non-Newtonian fluids cannot be predicted using single 
constitutive equation (unlike the case of viscous fluids). Therefore 
many models of non-Newtonian fluids are based either on 
"natural" modifications of established macroscopic theories or 
molecular considerations. The additional rheological parameter in 
the constitutive equations of non-Newtonian fluids are the main 
culprit for the lack of analytical solutions. The resulting equations 
are more complex and higher order than the Navier-Stokes 
equations. Hence these equations have been attracted from 
modelling as well as solutions point of view. The advancement 
in the study of non-Newtonian fluids has been made even by the 
recent investigators (See [1-10] and many studies therein). 

The flow induced over a stretching surface is very well 
documented problem in fluid mechanics. It is encountered in 
extrusion of polymer sheet from a die, glass fiber and paper 
production, continuous casting, cooling of metallic plate in a bath 



etc. Such flow problem in presence of heat transfer has been 
attracted by the researchers due to its applications in polymer 
processing technology. The quality of end product in industry 
depends upon both the stretching and cooling rates. Further, the 
thermal radiation effect has pivotal role in nuclear plants, gas 
turbines and devices for satellites, space vehicles, aircraft etc. The 
literature on this topic is quite sizeable. Rana and Bhargava [11] 
presented the numerical analysis for heat transfer of nanofluid over 
a nonlinearly stretching sheet. Bhattacharyya et al. [12] analyzed 
the solutions of boundary layer flow of viscoelastic fluid and heat 
transfer over a stretching sheet with internal heat generation or 
absorption. Makinde and Aziz [13] numerically studied the 
boundary layer flow of viscous nanofluid bounded by a stretching 
sheet. They considered the transport equation which includes the 
effects of Brownian motion and thermophoresis. Mandal and 
Mukhopadhyay [14] considered the boundary layer flow and heat 
transfer towards an exponentially stretching porous sheet embed- 
ded in a porous medium with variable surface heat flux. They 
found that the momentum and thermal boundary layer thickness 
decrease with increasing exponential parameter. Hayat et al. [15] 
examined the heat transfer in flow of second grade fluid over a 
stretching sheet. Thermal radiation effect in the boundary layer 
flow by stretching surface has been explored by Sajid and Hayat 
[16]. Bhattacharyya [17] discussed the unsteady stagnation point 
flow towards a stretching surface. Effect of heat transfer in flow 
over an exponentially stretching surface has been explored by 
Mukhopadhyay [18]. The radiation effect in flow of micropolar 
fluid towards a stretching surface is addressed by Hussain et al. 
[19]. Rashidi et al. [20] developed approximate solutions for heat 
transfer analysis in flow of micropolar fluid. Moreover, the interest 
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Figure 1. fi-curves of the functions f'(0) and 0 (0) at 10 th order 
of approximation. 
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in the study of magnetohydrodynamic flow for an electrically 
conducting fluid over heated surface is motivated by its great value 
in a wide range of engineering problems such as plasma studies, 
petroleum industries, MHD power generators, cooling of nuclear 
reactors, the boundary layer control in aerodynamics and crystal 
growth. Hence Turkyilmazoglu [21] found exact solution for 
magnetohydrodynamic flow of viscous fluid due to a rotating disk. 
Hayat and Nawaz [22] has investigated the Soret and Dufour 
effects in mixed convection three dimensional boundary flow of an 
electrically conducting second grade fluid over a vertical stretching 
sheet. Ahmad and Nazar [23] considered the problem of unsteady 
magnetohydrodynamic viscoelastic fluid flowing towards a stag- 
nation point on a vertical surface. Pal and Mondal [24] discussed 
the hydromagnetic flow of viscous fluid over a stretching surface in 
presence of both electric and magnetic fields. Abel et al. [25] 
presented MHD flow analysis for viscoelastic fluid. Both viscous 
and Ohmic dissipations are presented in this attempt. More, the 
analysis here is made when magnetic and electric fields are 
present. Hayat and Qasim [26] considered radiation effect in 
MHD flow of second grade fluid over unsteady porous stretching 
surface. The effect of internal heat generation in hydromagnetic 
non-Darcy flow and heat transfer over a stretching surface with 
thermal radiation and Ohmic dissipation is examined by 
Olanrewaju [27]. Elbashbeshy et al. [28] numerically analyzed 
the problem of unsteady laminar two-dimensional MHD bound- 
ary layer flow and heat transfer of an incompressible viscous fluid 

Table 1. Convergence of homotopy solutions when 



ai=0.02, £ = 0.05, a2 = 0.02, S = 
Re = 0.7, 7?i=0.3, M = 0.1, Pr = 


0. 2, A 

1, Ec = 


= 0.5, =0.3, 
= 0.5. 




Order of 

approximation — f"(0) 




-0'(O) 


1 1.0419 




1 .0059 


2 1.0720 




1 .0079 


5 1.1210 




1 .0041 


10 1.1442 




0.99450 


12 1.1458 




0.99211 


14 1.1458 




0.99051 


40 1.1458 




0.99051 


doi:10.1371/journal.pone.0083153.t001 
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Figure 2. Influence of M on f'(r\). 
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over a porous surface in the presence of thermal radiation and 
internal heat generation or absorption. MHD flow caused by a 
rotating disk is presented by Rashidi et al. [29]. The well-known 
Jeffery-Hamel problem in presence of magnetic field is examined 
by Motsa et al. [30] . Most of the studies on MHD flow over a 
stretching surface with heat transfer do not take into account the 
effect of electric field and Ohmic dissipation. Very little exists yet 
about such aspects in the stretched flows of viscous fluids. Such 
consideration further narrowed down when non-Newtonian fluids 
have been considered. To our knowledge there is only one such 
attempt for viscoelastic fluid [25]. The fluid employed although 
exhibits the normal stress effects but it cannot describe the features 
of shear thinning or shear thickening. Having such in view, the 
flow of third grade fluid is considered. This fluid even can capture 
shear thinning/shear thickening effects for one-dimensional flow 
over a rigid surface. The main objective here is to analyze the two- 
dimensional flow of third grade fluid over an unsteady stretching 
sheet. The effects of both electric and magnetic fields are retained 
in the momentum and energy equations. Thermal radiation and 
Ohmic dissipation are taken into account. The solutions for 
velocity component and temperature are developed by homotopy 
analysis method (HAM) [31-40]. The plots of physical quantities 
of interest reflecting the novel features of embedded parameters in 
the problems are given and analyzed. Tables for skin friction 
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Figure 3. Influence of p on f'Cn). 
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Figure 4. Influence of on f'(r\). 
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Figure 6. Influence of A on f'(r\). 
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coefficient and local Nusselt number are made and explained 
carefully. 

Mathematical Formulation 

We examine the two-dimensional boundary layer flow of 
magnetohydrodynamic (MHD) third grade fluid over a porous 
stretching surface. Here the fluid is electrically conducting in the 
presence of applied magnetic B = (0, Bq,0) and electric 
E= (0, 0, — Eq) fields. The flow is because of stretching of sheet 
from a slit through two equal and opposite forces. The sheet 
velocity is taken linear parallel to the flow direction. The electric 

and magnetic fields obey the Ohm's law J = o(^E-\- V x B^j. Here 

J is the Joule current, o is the electrical conductivity and V is the 
fluid velocity. The induced magnetic field and Hall current effects 
are ignored subject to small magnetic Reynolds number. Both the 
electric and magnetic fields contribute into the momentum and 
thermal boundary layer equations. The relevant equations in the 
aforestated conditions can be expressed as follows: 



du dv q 

dx dy 
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du 



d u 



EqBq — B 0 u) 



fdT dT 
d 2 T (du 



dT 



du d 2 u du d 2 u dud 2 u 
dy dtdy dy dxdy dy dy 2 



(3) 



In above equations u and v denote the velocity components in the 
x and y directions, a*, and are the fluid parameters, v is the 
kinematic viscosity, p is the density of fluid, T is the fluid 
temperature, K is the thermal conductivity of fluid, c p is the 



E f =/W =0.2, p=cci=0.1, A=1,S=0.4.Re=0.5 



Re=A=a2=0.2,ai=M=0.1,Ei=p=0.5 





Figure 5. Influence of a 2 on f'(r\). 
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Figure 7. Influence of S on f'(r|). 
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Figure 10. Influence of M on 0(r|). 
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specific heat at constant pressure and the radiative heat flux q r is 
first given by Sparrow and Cess [41] and Raptis [42] 



4a* dT 4 



(4) 



3*i dy ' 

where cr* is the Stefan-Boltzmann constant and k\ is the mean 
absorption coefficient. Through expansion of r 4 ^47^r-37^, 
Eq. (3) becomes 



(dT t ^dT dT 

dx dy 



pCp Vdi JrU ~ JrV 



3/c, 



d L T 



+ ai 



du d 2 u du d 2 u dud 2 u 



-2j8 3 



dy dtdy dy dxdy dy dy 2 
+ (j(uB 0 -E 0 ) 2 . 

The subjected conditions can be mentioned as follows: 
u(x, 0) = U w , v(x, 0) = V w , T(x, 0) = T w , 



(5) 



m->0, T^>T W , as j->oo, 
with F B . defined by 

Vo 



(l-ct) 



1/2' 



(6) 



(7) 



Here the mass transfer at surface with V w <0 is for injection and 
V w >0 for suction. Also the stretching velocity U w (x,t) and the 
surface temperature T w (x, t) are taken in the forms: 



(XX ax 
U w (x,t)=j-^ 9 T w (x,t) = T o0 + T Q — ^, (8) 



2v(l-cf) 



where a and c are the constants with a>0 and c>0 (i.e ct< 1). 
If i/f is the stream function then defining 



i=J^;y^=V^m,o=^^, (9) 
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Figure 9. Influence of E n on f'[r\). 
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Figure 11. Influence of p on 0(r|). 
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Figure 12. Influence of on 0(r|). 
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Figure 14. Influence of A on 6(tj). 
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The incompressibility condition is identically satisfied and the 
resulting problems for / and 6 are reduced into the following forms 
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(14) 



Here Re denotes the Reynolds number, M the magnetic 
parameter, E\ is the electric parameter, oc\ and 0C2 and are the 
fluid parameters, A is the suction parameter, S is the unsteadiness 
parameter, is the radiation parameter, Pr is the Prandtl 
number and Ec is the Eckert number. 
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The local skin friction coefficient is defined as 
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The velocity and temperature can be expressed in the set of base 
functions 
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where a m n and b m n are the coefficients. 

The initial guesses fo and 6o in homotopy solutions are taken 
through the expressions 




Figure 17. Influence of R d on Q(r\). 
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The auxiliary linear operators and their associated properties are 



im d ' f df i(ff\ d2 ° 



(21) 



satisfy the following properties 

L/[Ci + C 2 expfa) + C 3 exp( - rj)} = 0, (22) 



L 0 [C 4 exp(^) + C 5 exp( - ^)] = 0, 

Where C/(/ = 1 — 5) depict the arbitrary constants. 
The zeroth order problems are 



(23) 



(i - p )L f \f(^p)-Mr,)] = P h f n f [f{n,P)] , 

(1 -/OM§fa,/O-0ofa)] =/^N fl [jW), 

a 1 =a 2 =Rd=M=0.2,p=0.1,A=1,S=0.4,Pr=1,Re=Ec=0.5 



(24) 



a* 



a 4 




Ef = o.o 

E 1 = 0.2 
E 1 = 0.4 
E 1 = 0.6 
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Figure 20. Influence of Ec on 0(t|). 
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Table 2. Numerical values of skin friction coefficients Re l J 2 Cf for different values of physical parameters. 
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The general solutions of the Eqs. (31)— (32) are 

fmin) =fmin) + Ci + C 2 exp(*/) + C 3 exp( - rj), 
O m (rj) = 0* m (rj) + C 4 exp(^) + C 5 exp( - rj), 

in which f^iyj) and Q* m (?]) denote the special solutions. 

Convergence of the Derived Solutions 

We note that the series solutions (33) and (34) contain the non- 
zero auxiliary parameters hf and he. These parameters are useful 
in adjusting and controlling the convergence. The hf and 
ho — curves are plotted for 1 0 1 order of approximation in Fig. 1 
for the suitable ranges of the auxiliary parameters. Here the 
suitable values for hf and ho are — 1.5 <hf < —0.53, 
— 1.35 <ho < —0.4. Furthermore, convergence of series solution 
is checked and shown in Table 1. Note that the series solutions 
converge at 26th order of approximation up to 6 decimal places. 

Results and Discussion 

This section illustrates the impact of physical parameters. The 
results are displayed graphically in the Figs. 2-20. The conclusions 
for flow field and other physical quantities of interest are drawn. 
The numerical values of the skin friction coefficient and local 
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Table 3. Numerical values of Nusselt number Re T x / 2 Nu 


x for different values of physical parameters. 
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Nusselt number are presented in the Tables 2 and 3 for various 
values of o<2, S, M, E\, Re, Rj, Pr and Ec. Fig. 2 displays 
the effect of Hartman number M on velocity profile by keeping 
other physical parameter fixed. It is of interest to note that the 
velocity profile decreases with an increase in M whereas the 
boundary layer thickness reduces. Clearly by increasing magnetic 
force, the Lorentz force increases which cause resistance in the 
fluid flow and consequently the velocity profile decreases. Fig. 3 
shows the influence of third grade parameter on the velocity 
profile f'(jl). Here we noticed that the velocity increases near the 
wall with an increased whereas it vanishes away from the wall. 
Figs. 4 and 5 illustrate the variation of second grade parameters CL\ 
and 0C2 on the velocity profile f'(r\) respectively. It is observed that 
the velocity profile f'(f]) is an increasing function of 0L\. The 
velocity profile also increases when CL2 is increased. Fig. 6 is plotted 
for the effects of the suction parameter A on the velocity profile 
f f (fj). The velocity profile decreases by increasing parameter A 
and further the boundary layer is also decreasing function of A. 
Fig. 7 is sketched for the influence of unsteadiness parameter S on 
the velocity profile. The velocity profile and the thermal boundary 
layer decreases for larger values of S. The behavior of Reynolds 



number Re on velocity profile is shown in Fig. 8. It is observed 
that the velocity profile decreases with an increase in Reynold 
number Re. The influence of electric parameter E\ is shown in 
Fig. 9. This Fig explains that as the electric parameter E\ 
increases, the velocity boundary layer increases near the plate with 
small rate but increases away from the stretching plate more 
rapidly. In fact the Lorentz force (arising due to the electric field 
acts like an accelerating force) reduces the frictional resistance 
which causes to shift the stream line away from the stretching 
sheet. Fig. 1 0 portrays the effects of magnetic parameter M on the 
temperature profile 6(rj). It is depicted that temperature profile 
and thermal boundary layer thickness increase with an increase in 
magnetic parameter. Fig. 1 1 is the plot of temperature profile 6(fj) 
for various values of third grade parameter The effect of third 
grade parameter ft on 6(rj) shows a decrease near the wall. The 
boundary layer thickness also decreases. Figs. 12 and 13 describe 
the effects of second grade parameters ot\ and 0L2 on temperature 
profile 9(rj). Fig. 12. depicts that the effect of second grade 
parameter is to reduce the temperature distribution in the 
boundary layer which results in thinning of the boundary layer 
thickness. Same behavior is shown in Fig. 13 for various values of 
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Figure 21. Variation of velocity f'(r\) and shear stress f'(r|) with r\ for several values of Hartman number M. 
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0(2. The influence of suction parameter A and unsteadiness 
parameter S are analyzed in the Figs. 14 and 15. Here the 
temperature profile decreases with the increase of unsteadiness 
parameter S and the suction parameter A. Further the thermal 
boundary layer also decreases by increasing both the unsteadiness 
parameter S and the suction parameter A. Fig. 16 shows that the 
temperature profile and thermal boundary layer is decreasing 



function of Reynold number Re. The effects of thermal radiation 
parameter on temperature is shown in Fig. 17. It is revealed 
that the radiation parameter Rj causes increase in the fluid 
temperature 0(ff). On the other hand the thermal boundary layer 
thickness also increases. In Fig. 18 the influence of electric 
parameter E\ on temperature profile is given. This Fig. depicts 
that the temperature profile and the boundary layer thickness 
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Figure 22. Variation of velocity f (r|) and shear stress f'(r|) with x\ for several values of electric parameter E v 
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increase with an increase of electric parameter E\. Fig. 19. 
illustrates the effects of Prandtl number Pr on the temperature 
profile 6(rj). Both the temperature and thermal boundary layer 
thickness are decreased by increasing Pr. We displayed the 
temperature field for various values of Eckert number Ec in 
Fig. 20. The effect of Eckert number is to increase the temperature 
boundary layer thickness due to the frictional heating. Fig. 21 
shows the effects of Hartman number M on velocity f'(rj) and 
shear stress f"(r\). With the increase in M, the velocity field f'(r\) 
decreases near the wall and vanishes far away from the wall while 
shear stress f"(rj) has same behavior for larger values of Hartman 
number M. An opposite behavior is noted when 0 <f] <0.6. 
Fig. 22 demonstrates the effects of electric parameter E\ on 
velocity f'(j\) and shear stress f"(tj). It is worthmentioning to 
point out that velocity is increasing function of electric parameter 
E\ near the wall whereas opposite behavior for shear stress is 
observed for 0 < f] < 1 . The numerical values of skin friction 
coefficient for various physical parameters are shown in Table 2. 
Here the magnitude of skin friction coefficient increases with the 
increase of second grade parameters (tx\, 0(2), third grade 
parameter unsteadiness parameter S, Hartman number M 
and Reynold number Re whereas it decreases with an increase in 
electric parameter E\. Table 3 shows the effect of physical 
parameters on heat transfer characteristics at the wall —0'(Q). 
From this table, we observe that for large values of second grade 
parameters (<x\, 0C2), third grade parameter unsteadiness 
parameter S, radiation parameter Rj and Prandtl number Pr 
the heat transfer coefficient at the wall — 0'(O) increases while it 
decreases for Hartman number M, Reynold number Re, electric 
parameter Ej and Eckert number Ec. 
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